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Abstract. The main goal of this paper is to study properties of the iterated 
integrals of modular forms in the upper halfplane, eventually multiplied by z s ~ l , 
along geodesies connecting two cusps. This setting generalizes simultaneously the 
theory of modular symbols and that of multiple zeta values. 



§0. Introduction and summary 

This paper was inspired by two sources: theory of multiple zeta values on the 
one hand (see [Za2]), and theory of modular symbols and periods of cusp forms, 
on the other ([Mai], [Ma2], [Shl]-[Sh3], [Me]). Roughly speaking, it extends the 
theory of periods of modular forms replacing integration along geodesies in the 
upper complex half-plane by iterated integration. Here are some details. 

0.1. Multiple zeta values. They are the numbers given by the /c-multiple 
Dirichlet series 

C(m 1 ,...,m k ) = V — ^ (0.1) 

' . . .rii 



. b i . . . / 1 1 

0<ni<---<n fc 1 k 



which converge for all integer rrii > 1 and nik > 1, or equivalently by the m-multiple 
iterated integrals, m = mi + ■ ■ ■ + m^, 

1 dzi f Zl dz 2 f Z2 [ Zm k-i dz mk 



C(m 1 ,...,m k )= ==i / =lf ... (0.2) 

JO z l JO z 2 Jo Jo 1 — z m k 

where the sequence of differential forms in the iterated integral consists of consec- 
utive subsequences of the form — , . . . , — , of lengths mfc, mk-i, • • • , rn\. 

z z 1 — z 

Easy combinatorial considerations allow one to express in two different ways 
products C(h, . . . ■ C(^i, • • • , irik) as linear combinations of multiple zeta values. 

If one uses for this the integral representation (0.2), one gets a sum over shuf- 
fles which enumerate the simplices of highest dimension occurring in the natural 
simplicial decomposition of the product of two integration simplices. 
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If one uses instead (0.1), one gets sums over shuffles with repetitions which 
enumerate some simplices of lower dimension as well. 

These relations and their consequences are called double shuffle relations. Both 
types of relations can be succinctly written down in terms of formal series on free 
noncommutating generators. One can include in these relations regularized multiple 
zeta values for arguments where the convergence of (0.1), (0.2) fails. 

For a very clear and systematic exposition of these results, see [De] and [Ral], 
[Ra2]. 

In fact, the formal generating series for (regularized) iterated integrals (0.2) ap- 
peared in the famous Drinfeld paper [Dr2], essentially as the Drinfeld associator, 
and more relations for multiple zeta values were implicitly deduced there. The 
question about interdependence of (double) shuffle and associator relations does 
not seem to be settled at the moment of writing this: cf. [Ra3]. The problem of 
completeness of these systems of relations is equivalent to some difficult transcen- 
dence questions. 

Multiple zeta values are interesting, because they and their generalizations ap- 
pear in many different contexts involving mixed Tate motives ([DeGo], [T]), defor- 
mation quantization ([Kon]), knot invariants etc. 

0.2. Modular symbols and periods of modular forms. Let T be a congru- 
ence subgroup of the modular group acting upon the union H of the upper complex 
half-plane H and the set of cusps P 1 (Q). 

The quotient Y \ H is the modular curve X-p- Differentials of the first kind on 
Xr lift to the cusp forms of weight 2 on H (multiplied by dz). 

The modular symbols {a,/3}r G Hi(Xr,Q), where a, (3 G P 1 (Q), were intro- 
duced in [Mai] as linear functionals on the space of differentials of the first kind 
obtained by lifting and integrating. The fact that one lands in Hi(Xr, Q) and not 
just Hi(Xr,~R) is not obvious. It was proved in [Drl] by refining a weaker argu- 
ment given in [Mai] . This is equivalent to the statement that difference of any two 
cusps in Y has finite order in the Jacobian, or else that the mixed Hodge structure 
on H 1 (Xr \ {cusps}, Q) is split (cf. [El]). One of the basic new insights of [Mai] 
consisted in the realization that studying the action of Hecke operators on modular 
symbols one gets new arithmetic facts about periods and Fourier coefficients of cusp 
forms of weight two. 

The further generalizations of modular symbols proceeded, in particular, in the 
following directions. 
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(a) In [Ma2] it was demonstrated that the same technique applies to the integrals 
of cusp forms of higher weight, eventually multiplied by polynomials in z, producing 
the similar information about their periods and Fourier coefficients. In principle, 
such integrals cannot be pushed down to X-p, but they can be pushed down to 
the appropriate Kuga-Sato varieties over Xp that is, relative Cartesian powers of 
the universal elliptic curve. In this way, modular symbols of higher weight can be 
interpreted as rational homology classes of middle dimension of Kuga-Sato varieties: 
cf. [Shl]-[Sh3]. 

(b) Pushing down an oriented geodesic connecting two cusps in H to Xp, we 
get a singular chain with a boundary in cusps, which is a relative cycle modulo 
cusps with integral coefficients. This is the viewpoint of [Me]. Hence it is more 
natural to consider the relative/non-compact version of modular symbols, and allow 
integration of the Eisenstein series, that is, differential forms of the third kind with 
poles at cusps as well. The same remark applies to the modular symbols of higher 
weight. 

This refinement appears as well in the study of the "noncommutative boundary" 
of the modular space, that is, the (tower of) space(s) T \ P 1 (R), cf. [MaMar]. 
Namely, it turns out that the relative 1-homology modulo cusps (and additional 
groups of similar nature) can be interpreted as (sub) groups of the K-theory of the 
noncommutative boundary. 

In this paper I suggest a generalization in the third direction, namely 

(c) The study of iterated integrals of cusp forms and Eisenstein series, eventually 
multiplied by a power of z, along geodesies connecting two cusps. Some of these 
integrals can be pushed down to Xp and thus produce a de Rham version of modular 
symbols which assigns iterated (eventually regularized) periods to the elements of 
the fundamental groupoid of (Xp, {cusps}) instead of its 1-homology group. One 
may call them noncommutative modular symbols. 

Other integrals can only be pushed down to the Kuga-Sato varieties, or prefer- 
ably, to some (covers of the) moduli spaces Mi >n , in the same vein as it was done for 
multiple zeta values and M 0)n in [GoMa]. The related geometry deserves further 
study, both for integrands related to cusp forms and to Eisenstein series. 

Notice in conclusion that the discussion above implicitly referred only to the case 
of <SX 2 -modular symbols. It would be quite interesting to extend it to groups of 
higher rank, along the lines of [AB] and [AR]. 

0.3. Summary of this paper. I recall the basic properties of iterated integrals 
of holomorphic 1-forms on a simply connected Riemannian surface in §1. The shuf- 
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fle relations for the iterated integrals are reflected directly in terms of a generating 
function J stating that it is a group-like element with respect to a comultiplication, 
cf. Proposition 1.3.1. 

Then I turn to the main object of study. In §2 I define 1-forms of modular and 
cusp modular type, introduce and study the iterated and total Mellin transform for 
families of such forms. The functional equation for the total Mellin transform is 
deduced which extends the classical functional equation for L-series. 

Using only critical values of these Mellin transforms, I introduce in 2.6 an iterated 
modular symbol as a certain noncommutative 1-cohomology class of the relevant 
subgroup of the modular group. 

In §3, I study the representation of such Mellin transforms at integer values of 
their Mellin arguments in terms of multiple Dirichlet series. The results differ from 
the classical ones expressed by the identity (0.1) = (0.2) in two essential respects. 
First, iterated integrals are only linear combinations of certain multiple Dirichlet 
series. Second, the latter are not of the usual type 



in fact, their coefficients depend on pairwise differences rij — n\. 

In §4, the properties of the multiple Dirichlet series which emerged in §3, are 
axiomatized, and the shuffle relations for them are deduced. This requires, however, 
a considerable extension of the initial supply of series; the system of those coming 
from 1-forms of modular type is not closed. 

§5 is dedicated to the iterated analogs of the so called Eichler-Shimura and 
Manin relations for periods of cusp forms. Whereas the relations of the first type are 
quite straightforward, the relations of the second type, involving Hecke operators, 
are not obvious. The results presented here (Theorem 5.3) are preliminary, they 
clearly allow generalizations and deserve further study. 

Finally, in §6 I return to the formalism of §1 and extend it by allowing our 
integrands to have logarithmic singularities at the boundary. A version of the 
regularization procedure I use here is the same as in Drinfeld's paper [Dr2]. It 
exploits complex analyticity in place of Boutet de Monvel's technique of [De] and 



Using the Manin-Drinfeld theorem on cusps, I suggest a generalization of Drin- 
feld's associator and extend to this case a part of the identities satisfied by the 



0<n 1 <---<n k 





[Ra2]. 
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latter. This list includes the group-like property, the duality, and the hexagonal 
relation, which turn out to have the same source as the Shimura - Eichler relations 
for the periods of cusp forms. To the contrary, the pentagonal relation seems to be 
specific for the original Drinfeld's associator. 

Acknowledgements. I am grateful to A. Levin and A. Goncharov who have read 
the first draft of this paper and made a number of useful suggestions which are 
incorporated in the text. 

§1. Iterated integrals of holomorphic 1— forms 

1.1. Setup. Let X be a connected Riemann surface, not necessarily compact, 
Ox its structure sheaf of holomorphic functions, VL X the sheaf of holomorphic 1- 
forms. If a; is a (local) 1-form, z G X a point, uj(z) denotes the value of uj at z, i.e. 
the respective cotangent vector. 

Let V be a finite set which will be used as a set indexing various families. 
Consider the completed unital semigroup ring freely generated by V. We will write 
it as the ring of associative formal series C{{Ay)) where Ay := (A v \ v G V) are 
noncommuting free formal variables. 

More generally, we may consider the ring Ox(U)((A v )) where Ox(U) is the ring 
of holomorphic functions on an open subset U C X (A v commute with Ox(U)), 
and the bimodule Q. x (U)((Ay)) over this ring, connected by the differential d such 
that dA v = for all v G V. Varying U, we will get two presheaves; the sheaves 
associated with these presheaves are denoted Ox {{Ay)), resp. Q x ((Ay)), and d 
extends to them, so that Ker<i is the constant sheaf C((Ay)}. 

Let uy := (uj v \ v G V) be a family of 1-forms holomorphic in U and indexed by 
V. Put 

O := ^ A v^v ■ (1.1) 
vev 

The total iterated integral of this form along a piecewise smooth path 7 : [0, 1] — > U 
is denoted J 7 (f2) or J 1 (uy) and is defined by the following formula: 

00 rti rtn-i 

j 7 (q):=i+]t/ rmti) i*(m2>-- / nmn) e c((a v )) (1.2) 

n=1 Jo Jo Jo 

where the integration is taken over the simplex < t n < • • • < t\ < 1. If 7, 7' with 
the same ends are homotopic, J 7 (0) = J 7 '(0). 
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Putting Zi = 7(tj) G X, a = 7(0), z = 7(1), and considering the whole integral 
as a function of a variable z we will also write (1.2) in the form 

00 fZ pZ\ fZn-l 

r a {Q) = r a (u v ) = i + J2 n (*0 / ■ ■ / «M • (i-3) 

If £7 is connected and simply connected, this expression is an unambiguously defined 
element of Ox(U)((Av)). Otherwise it is a multivalued function of z in this domain. 

The following result is classical. 

1.2. Proposition, (i) J„(fi) as a function of z satisfies the equation 

e^(Q)=fi(z)J a 2 (Q). (1.4) 

In other words, J„(fi) is a horizontal (multisection of the flat connection Vn := 
d — Iq on Ox{{Ay)) , where Iq is the operator of left multiplication by O. 

(ii) If U is a simply connected neighborhood of a, J„(0) is the only horizontal 
section with initial condition J% = 1. Any other horizontal section K z can be 
uniquely written in the form J Z (Q)C, C G C{{Ay)). In particular, for any b G U , 

jm = J z a {v)jm (1.5) 



Proof, (i) follows directly from (1.3). Since J z (£l) is an invertible element 
of the ring Ox(U)((Ay)) , we can form J Z (Q,)~ 1 K Z and then directly check that 
d{ J z (fl)~ 1 K z ) = 0. Hence this element belongs to C((Ay)), and moreover equals 
its value at z = a that is, K a . Choosing K z = J£(Q), we get (1.5). 

1.3. Jn(n) as a generating series. Clearly, we have 

00 

J z a (oj v ) = J z a (Q) = l + Y, A Vl ...A Vn I z (u; Vl ,...,uj Vn ), (1.6) 

n=l (v!,...,v n )&V n 

where 

rz rzi rz n -i 

I z (u Vl ,...,uj Vn ) = u Vl (z\) I u V2 (z 2 )--- uj Vn (z n ) (1.7) 

J a J a J a 



are the usual iterated integrals. 
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In the remaining part of this section, and in the main body of the paper, we 
will encode various (infinite families of) relations among the iterated integrals (1.7) 
in the form of relations between the generating functions J^(uv)- Generally, our 
relations between the generating functions will be (noncommutative) polynomial 
ones. They may also involve different families (coy), different integration paths, 
and some linear transformations of the formal variables A v : cf . especially Theorem 
2.3, Proposition 5.1.1, Theorem 5.3, and subsection 6.5 (in the context requiring a 
regularization) . 

1.4. Basic relations between total iterated integrals. There are three 
types of basic relations, which we will call group-like property, cyclicity, and func- 
toriality respectively. 

1.4.1. Proposition. Consider the comultiplication 

A : C((A V )) -> C((A V ))® C C((A V )), A(A V ) = A v <g> 1 + 1 <g> A v 
and extend it to the series with coefficients Ox and Q x . Then 

A (J z a (cj v )) = J z a {uj v ) ® 0x J z a (uv) ■ (1.8) 

Proof. Both sides of (1.8) satisfy the equation dJ = A(f2) J and have the initial 
value 1 at z = a. 

NB. Coefficientwise, (1.8) is a compact version of shuffle relations for iterated 
integrals (1.7). 

1.4.2. Cyclicity. Let 7 be a closed oriented contractible contour in U , ai, . . . , a n 
points along this contour (cyclically) ordered compatibly with orientation. Then 

J£(n)J%(n)...JZ-*(ci)JZ(ci) = 1. (1.9) 

This follows from (1.5) by induction. 

1.4.3. Functoriality. Consider an automorphism g : X — > X such that g* 
maps into itself the linear space spanned by ui v . In particular, there is a con- 
stant matrix G = (g vu ) with rows and columns labeled by V such that g*(w v ) = 
^2 u 9vu^u- Define the automorphism g* of any of the ring/module of formal se- 
ries C{{A V )), C(X){{A V )), 1 (X)((A y )) by the formula g*{A u ) = J2 v A v g vu . On 
coefficients g* acts identically. 
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1.4.4. Claim. We have 

J£lw) = g.{J* a {uv)) . (1.10) 

Proof. In fact, both sides coincide with J*(g*(uv))- We will give below a 
calculation which proves a slightly more general statement. 

1.5. A variant: multiple lower integration limits. Somewhat more gener- 
ally, in the simply connected case we can consider a family of points (a,) := (a^„) 
in A indexed by pairs i = 1,2, 3, v e V. 

Given such a family and ujy, we can construct the following formal series in 
C(X){{Ay)) with constant term 1: 

oo 

n=o (tii,...,t)„)ey» 

where 

I Zi, vl ,...,a n , Vn (f*>v 1 ,---,Wv n )--= / U V2 (Z 2 )--- U Vn (z n ). 

(1.12) 

As above, 2; G A denotes a variable point, the argument of our functions. Then we 
have 

^fa.)(^)=^(a.)(^V) (1.13) 

and 

J(a.)(!Ov) = J Z a("v) J^V) ■ (1.14) 

The series (1.11) satisfies the following functoriality relation generalizing (1.10): 
1.5.1. Claim. We have 

Proof. We will check that both sides coincide with J* a ^(g*(uiy)). In fact, 
JJJ v(z) = z/((7^) so that, removing g step by step from the integration lim- 
its, we get 

i |ai,„ 1 ,... )S a n ,„„ ( w «i» • • • > w O = I a 1 , vl ,...,a n , Vn (0*( w «i)» • • • » 0*( W «J) • 
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Multiplying the l.h.s. by A Vl . . . A Vn and summing, we get the l.h.s. of (1.15). 
On the other hand, 

Yl A ^ ■ ■■ A vJa 1 , vl ,...,a n , Vn (9*(Uv 1 ), ■ ..,g*(w Vn )) = 
V!,...,v n €V n 

E ^ v:,„ *„,„,( E *.*<■«. E owo = 

v 1 ,...,v n ev n mev u n ev 

= A Vl g VljUl . . . A Vn g VnjUn I ai v ^ ,...,a n , Vn Kn • • • > = 

wi,...,u„ev™ 
ui,...,u n ev" 

= £* ^ ^■••V« 1 ,„ 1 ,...,a„, % Kv,W»J 

\»i,-,v n 6V" 

Summation over n produces the r.h.s. of (1.15), proving the lemma. 

1.6. A variant: nonlinear O. Let now O G 0^(£/") ((^4y)) be an arbitrary 
form without a constant term in A„: 

oo 

V = J2 E ^•••^A,...,^, (1.16) 

n=l («i,...,v„)€V« 

where ft^, ...,„„ G ^(t/). 

The total iterated integrals J 7 (0) and J„(0) are defined by exactly the same 
formulas (1.2) and (1.3). It is not true anymore that the coefficients of this series 
are the usual iterated integrals. However, an analog of the Proposition 1.2 and the 
cyclic identity remain true: 

1.6.1. Proposition. J„(0) as a function of z satisfies the equation 

dj*(n) = n(z)j*(n). (i.i7) 

If U is a simply connected neighborhood of a, J„(0) is the only horizontal section 
with initial condition J% = 1. Any other horizontal section K z can be uniquely 
written in the form J*(Q)C, C G C{{Ay))- In particular, for any b G U , 

jm = J z a {v)jm (l.is) 
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1.6.2. Corollary. Let 7 be a closed oriented contractible contour inU , a±, . . . , a n 
points along this contour (cyclically) ordered compatibly with orientation. Then 

J2(fl)J2(fl)...^-'(n)jj(n) = 1. (1.19) 



Notice in conclusion that the integral formula (0.2) for the multiple zeta values is 
not quite covered by the formalism reviewed so far because the integrands in (0.2) 
have logarithmic poles at the boundary. We will return to this situation in §6, to 
which some readers may prefer to turn right away. However, for applications to the 
integration of cusp forms in §2 - §5 the regular case treated here suffices. 

§2. 1— forms of modular type, iterated Mellin transform, 
and noncommutative modular symbols 

2.1. Setup. In this section, X will be the upper half plane H and z the 
standard complex coordinate. H is endowed with the metric of constant curvature 
— 1: ds 2 = \dz\ 2 1 '(Im z) 2 . 

The limits of integration in our iterated integrals generally lie in H, but may be 
"improper" as well, that is, belong to the set of cusps Q U {zoo}. If this is the case, 
we always assume that the respective integration path in some neighborhood of the 
cusp coincides with a segment of a geodesic curve. 

Our 1-forms generally will have the following structure. 

2.1.1. Definition, (i) A 1-form u on H is called a form of modular type, if 
it can be represented as f(z)z s ~ 1 dz where s is a complex number, and f(z) is a 
modular form of some weight with respect to a congruence subgroup of the modular 
group. 

The modular form f(z) is then well defined and called the associated modular 
form (to ui), and the number s is called the Mellin argument of u. 

(ii) u> is called a form of cusp modular type if the associated f(z) is a cusp form. 

To fix notation, we will recall below some classical facts. 

2.1.2. Action of automorphisms. Any matrix 7 G GL^(R) defines a holo- 
morphic isometry of if, namely zh [7] 2 where [7] is the fractional linear transfor- 
mation corresponding to 7. We will denote this automorphism also 7. It induces 
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the inverse image maps on the sheaves (0^) (glr of holomorphic tensor differentials 
of degree r: 

7 *(/(z) (dzY) = /([ 7 ]z) (d[ 7 ]z) r = (det 7 )7([7]^) (J^) 2 - ^ 

where (c 7 , d 7 ) is the lower row of 7. 

If one identifies (fi^)® r with by sending (dz) r to 1, (2.1) turns into the 
action of weight 2r on functions which is traditionally written as a right action: 

/|[ 7 ]2r(*) := (det 7 )V([7]^)(c 7 ^ + rf 7 )- 2r . (2.2) 

Assume that f(z) (dz) r is invariant with respect to 7 . Then, writing f(z)z s ~ 1 dz = 
f(z)(dz) r ■ z s ~ 1 (dz) 1 ~ r , we see that 

7*(/(z) z fl_1 dz) = (det 7 ) 1 "7(z)(a 7 2 + 6 7 ) s " 1 (c 7 2 + d^'^'dz. (2.3) 

where (a 7 , 6 7 ) is the upper row of 7 . In particular, if 2r > 2 is an integer, 7 * maps 
into itself the space of 1-forms spanned by 

f(z)z s - 1 dz, 1 < s < 2r-l, s e Z. (2.4) 

More generally, if 

7 *(/(z)(cte) r )=x(7)/WW r (2-5) 

for some x( 7 ) G C, then 

7 *(/(^)^- 1 ^) = (det 7 ) 1 " r f(z) X (l)(a 1 z + b i y- 1 (c 1 z + d 1 ) 2r - 1 - s dz 1 (2.6) 

and the space (2.4) will still remain invariant. 

We can apply this formalism to the spaces of modular forms of weight 2r 
with respect to a congruence subgroup V of (Z), i.e. to the functions / in 
(dz)~ r ((fi^)® r ) r . Two special cases will be of particular interest: 

(i) For any such /, the space of 1-forms spanned by (2.4) is T-invariant. 

(ii) Assume that F — F (N). This group is normalized by the involution 
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Therefore this involution maps into itself the space of r (iV)-modular forms, and 
the latter has a basis consisting of forms with 

g* N (f(z)(dzy) = e f f(z)(dz) r , e f = ±1. (2.8) 

Applying (2.6) with 7 = g^ we get, for any complex s, 

g* N (f(z)z s - 1 dz) = e f N r - s f(z)z 2r - 1 - s dz. (2.9) 

2.1.3. Geodesies and cusp forms. The geodesic from to too is the upper 
half of the pure imaginary line. The unoriented distance of a point iy on it to z is 
I log 2/ 1 - The exponential of this distance is thus y, if y > 1, and y~ , if y < 1. If we 
replace i by another reference point, even outside of imaginary axis, the exponential 
of the distance will behave like e°^y (resp. e ^ 1 )?/ -1 ) as y — > 00 (resp. y — > 0.) 

Let /(z) be a cusp form of weight 2r for a congruence subgroup. Then it can be 
represented by a Fourier series f(z) = Xl^Li c n e 2lTiriZ / N for some iV e Z+, whose 
coefficients are polynomially bounded: c n = 0(n c ) for some C > 0. Therefore we 
have \f(iy)\ = 0(e~ ay ) for some a > as y — > 00. From the previous analysis 
it follows that more generally, for any cusp form and any geodesic connecting two 
cusps, \f(z)\ = 0{e~ ay ^) for some a > as z tends along the geodesies to one of 
its ends, where this time y(z) means the exponentiated geodesic distance from z to 
any reference point in H, fixed once and for all. 

Let now uj(z) = f(z)z s ~ 1 dz be an 1-form of cusp modular type. Then the 
estimates above show that the following expected properties indeed hold. 

a) As zq — > zoo along the imaginary axis, the family to of holomorphic func- 
tions of z in any bounded domain H converges absolutely and uniformly to a holo- 
morphic function of z which is denoted j ioo u. The same remains true, if one 
replaces zoo by 0. 

These integrals are holomorphic functions of the Mellin argument s of uj as well. 

b) The sum (J^ + J*°) oj does not depend on z in H and is denoted f ioo uj. As 
a function of s, it is called the classical Mellin transform of uj. 

Denote this classical transform A(/; s). Assume that / satisfies (2.8). Then we 
have the classical functional equation 



A(f;s) = -e f N r - s A(f;2r-s), 



(2.10) 
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because in view of (2.9) 

fO rioo rgN(0) 



/ u> = - / u = - / w = - / g* N (u;) = -e f N r ~ s / f(z); 

J ioo Jo Jq N (ioo) J ioo J ioo 



Another identity in the same vein uses the fact that i / y^N is the fixed point of gjy 
so that A(/; s) can be written as 

A(/;a)= / u- g* N {uj). (2.11) 

J ioo J ioo 

This allows one to use the Fourier expansions of f(z) and f\[gN]2r(z) in order 
to deduce series expansions for A(/; s) (notice that the Fourier expansions cannot 
be term-wise integrated near z = because the formal integration produces a 
divergent series). 

Now we can finally write down the analogs of these definitions and results for 
iterated integrals. 

2.2. Definition, (i) Let fi,. . . , fk be a finite sequence of cusp forms with respect 
to a congruence subgroup, u)j(z) := fj(z) z Sj ~ 1 dz. The iterated Mellin transform of 
(fj) is, by definition, 

M(f l ,...,f k ;s 1 ,...,s k ) := ^(wi, . . . , w k ) = 

pz\ rz n -i 

wi(zi) / uj 2 (z 2 )--- uj n (z n ) (2.12) 

too J ioo J ioo 

(ii) Let f v = (fv\v e V) be a finite family of cusp forms with respect to a con- 
gruence subgroup, sy = (s v \ v G V) a finite family of complex numbers, uy = 
where uj v (z) := f v (z) z Sv ~ 1 dz. The total Mellin transform of fy is, by definition, 

TM(fy-Sy) :=J^(u V ) = 

oo 

= A Vl . . . A Vn M(f Vl , . . . , f Vn ; s Vl , . . . , s Vn ) (2.13) 

n=0 (v 1 ,...,v n )€V n 

(cf. (L3)). 
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Below we will assume that the space spanned by all uj v is stable with respect to 
some g* N . Then as in 1.4.3 denote by G = (g vu ) the matrix of this action on (u> v ), 
and by g^* the action of the transposed matrix on the formal variables (A v ). 

For example, if (to v ) and (f v (dz) Tv ) respectively can be represented as a union of 
pairs of forms, corresponding to the left and right hand sides of (2.9), the matrix G 
consists of two by two antidiagonal blocks each of which after the classical Mellin 
transform produces a functional equation of the form (2.10). 

2.3. Theorem, (i) If the space spanned by all u v is stable with respect to some 
g* N , we have the following functional equation: 

J° 00 (^)=^(J° 00 (^))- 1 . (2.14) 

(ii) In the assumptions of the Definition 2.2 (ii), denote the weight of f v by 2r v 
and assume that f v is an eigenvector for g* N with eigenvalue e v . Then the total 
Mellin transform (2.13) satisfies 

TM(f v ; s v ) = g*(TM(f v ; 2r v - sy))- 1 (2.15) 

where g* multiplies each A v by e v N rv ~ Sv . 

Proof. This is a straightforward corollary of the definitions and formulas (1.9) 
and (1.10) as soon as one has checked that the latter formulas are applicable to the 
improper iterated integrals of the 1-forms of cusp modular type. 

This check is a routine matter, since at each step of an iterated integration we 
multiply the result of the previous step by a holomorphic function of the type 
f(z) z s ~ x which is bounded by 0(e~ ay ^) as in 2.1.3 above as z tends to or zoo. 

Notice in conclusion that no analog of the functional equation (2.11) can be 
written for the individual Mellin transforms (2.12), because applying g n to the 
integration limits in them we get an expression which is not a Mellin transform in 
our sense. Only putting them all together produces the necessary environment for 
replacing the overall minus sign at the r.h.s. of (2.10) by the overall exponent — 1 
at the r.h.s. of (2.15). 

A similar reasoning establishes the iterated analog of (2.11): 

2.4. Proposition. We have 

TM(f v ; s v ) = (gN.J^Mr'Jl^M ■ (2.16) 
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2.5. Pushing down iterated integrals. Let u be an 1-form of modular type 
whose associated modular form has weight 2 with respect to a subgroup Y of the 
modular group, and whose Mellin argument is 1. In this case uj is T-invariant so 
that it can be pushed down to an 1-form v on X^ := T\H. Instead of integrating 
u) along a path in H, we can integrate v along the push-down of this path to X^. 
If all uj v have this property, all relevant iterated integrals can be pushed down to 
X°. 

This argument admits a partial generalization to higher weights. Assume that 
the modular form associated with uj has weight 2r > 2, whereas its Mellin argument 
is an integer belonging to the critical strip (2.4), 1 < s < 2r — 1. In this case the 
relevant simple integral along, say, {200, 0} can be pushed down to the Kuga-Sato 

variety X^? r ^ which is the (2r — 2)-th fibered power of the universal elliptic curve 

over Xy or rather its compactified smooth model. However, on X r we obtain 
an integral of a holomorphic form uj of degree 2r — 1 over a relative cycle of the 
same dimension which is > 1. Therefore iterated "line" integrals of such forms on 
H cannot be directly translated into integrals of the same type on X^ r ~ 2 \ 

On the other hand, one can generally define Chen's iterated integrals of forms of 
arbitrary degree, say, uj v on X^ r 2 \ which take values in the space of differential 

forms on the path space PX^ r ~ 2 ^ and not just C: cf. reports [Ch] and [Ha], as well 
as references therein. Studying properties of such iterated integrals in the modular 
case like presents an interesting challenge. 

Here I will restrict myself to explaining how uj looks like and why its periods 
coincide with integrals of u> along geodesies. For more details, see [Shi] , [Sh2] , and 
especially [Sh3]. 

Denote by the semidirect product Tx (Z 2r ~ 2 xZ 2r ~ 2 ) acting upon HxC 2r ~ 2 
via 

(7; n,m) (z,() := ([n/]z; (c 7 z + d 7 ) _1 (C + zn + to)) . 

Here n = (m, . . . , n 2r -2), m = (mi, . . . , m 2r -2), C = (Ci> • • • 5 Qir-2), and nz = 
(n 1 z,...,n 2r -2z). 

If f(z) is a holomorphic modular form of weight 2r, then f(z)dzAd(iA- ■ -A<i(2r-2 
is a l?( r ) -invariant holomorphic volume form on H x C 2r_2 . Hence one can push it 
down to (a Zariski open subset of) the quotient \ (H x C 2r " 2 ) which is a Zariski 
open subset of the respective Kuga-Sato variety. Denote by uj the image of this 
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form. Notice that it is common for all 1-forms of modular type u> = f(z)z s x dz 
with different Mellin arguments s. 

A detailed analysis of singularities performed in [Sh2] , [Sh3] shows that the map 
f i — * uj induces an isomorphism of the space of cusp forms of weight 2r with the 
space of holomorphic volume forms on an appropriate smooth projective Kuga- 
Sato variety. (As I have already remarked in the Introduction, it would be useful 
to replace it by the base extension (M 1 ^ r -2)x r -) 

The dependence of the period of uj on the integration path and on the Mellin 
argument is reflected in the choice of the relative cycle over which we integrate cD. 

More precisely, let a, (3 G P 1 (Q) be two cusps in H and let p be a geodesic joining 
a to (3. Fix (rii) and (mi) as above. Construct a cubic singular cell p x (0, l) 2r_2 — ► 
H x C 2r_2 : (z, (ti)) i— > (z, (ti(zrii + mi))). Take the 5 , 2 r -2^symmetrization of this 
cell and push down the result to the Kuga-Sato variety. We will get a relative cycle 
whose homology class is Shokurov's higher modular symbol {a, (3; n,m}r- From 
this construction it is almost obvious that 



2r-2 

f(z) ^ (rijZ + mi) dz = / Q. 

i= l J {a,P;n,m} r 



The singular cube (0, l) 2r_2 may also be replaced by an evident singular simplex. 
This can be useful for transposing the results of [GoMa] to the genus one moduli 
spaces. 

2.6. Noncommutative modular symbols and continued fractions. I will 
define in this subsection a generalization of modular symbols involving iterated 
integrals and allowing a mixture of forms of different weights with respect to the 
same subgroup V of SL{2, Z). 

Let (coy) be a family of linearly independent 1-forms of cusp modular type whose 
Mellin arguments are integers lying in the respective critical strip as in (2.4). Let 
T be a subgroup of modular group acting on the space spanned by (uj v ) as in (2.3). 
Denote by IT the multiplicative group of power series in (A v ) with constant term 
1. Clearly, the map J i— > g*J (see 1.4.3.) defines the left action of V on IT. 

2.6.1. Proposition— Definition, (i) For each a G P 1 (Q), the map T — > II : 

7 l— *■ Jjaity i s a noncommutative 1-cocycle Q a in Z 1 (r,II). 

(ii) The cohomology class of ( a in if 1 (r,n) does not depend on the choice of a 
and is called the noncommutative modular symbol. 
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Proof. We have, omitting for brevity, and using (1.9), (1.10): 

ja _ ja jja _ ja t ja \ 

J ~fj3a — J -ya J ~ff3a ~ J ~fa !*\ J /3a) 

which means that ( a is an 1-cocycle. Moreover, if b is another cusp, 

J^a = JbJjbJ^a = Jb Jjb(l*(Jb )) 

that is, Ca and (b are homologous. 

Remark. Assume that the cusp forms associated with (u v ) span the sum of all 
spaces of cusp forms of certain weights, and for each weight and each cusp form, 
all admissible Mellin arguments actually occur. Then the linear in A v term of £a 
encodes all periods of the involved cusp forms along all classical modular symbols 
corresponding to loops in X-p(C) starting and ending at the cusp Fa. 

2.6.2. Iterated integrals between arbitrary cusps. The group T generally 
does not act transitively on cusps, so that the components of cocycles ( a do not 
contain iterated integrals along all geodesies connecting two cusps. One can use 
the technique of continued fractions as in [Mai] , [Ma2] in order to express all such 
integrals through a finite number of them. 

Namely, choose a set of representatives C of left cosets V \ SL 2 (Z). Call the 
iterated integrals of the form («^( l °2))) ±1 ' 9 e P r ^ m ^ ve ones. Notice that when 
g T the space spanned by (u v ) is not generally g*-stah\e so that we cannot define 
9*- 

2.6.3. Proposition. Each J£ can be expressed as a noncommutative monomial 
in 7*(JJ) where 7 runs over V and J% runs over primitive integrals. 

Proof. First, we can write = {J^°)~ l Jl°° ■ So it remains to find a required 
expression for Jf^. Assume that a > 0; the case a < can be treated similarly. 
Consider the consequent convergents to a: 

_ Pn Pn-l P0 _P0 P-l _ 1 

q n q n - X q 1 q_i 

Put 

9k -\q k (-l) fc "V-J' k ~°--- n - 
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We have = gk( a ) £ SL 2 (2, Z). Put gj~ = ^kCk where 7^ G Y and G C are two 
sequences of matrices depending on a. Then from cyclicity we get 

n 

ja _ TT jPk/qk 

J ioo - 11 J p k _ 1 /q k _ 1 
k=0 

and from functoriality we obtain 

jPk/Qk _ I TC fe (ioo)N 

J W-i/ft-i~' fe *^ci(0) 

§3. Values of iterated Mellin transforms at integer points 
and multiple Dirichlet series 

In this section, we collect some formulas expressing iterated Mellin transforms 
(2.12) at integer values of their Mellin arguments as linear combinations of "multiple 
Dirichlet series". 

3.1. Notation. Consider a family of 1-forms u)y , v e V, satisfying the following 
conditions. First, 

00 

u v {z) = Yl c v , n e 27rinz z m *- 1 dz, c v , n G C , m v G Z, m v > 1. (3.1) 

n=l 

Moreover, assume that c„ >n = (9(n c ) for some C and each v. 

Until a problem of analytic continuation arises, we do not have to assume mod- 
ularity. The notation m v replacing former s v is chosen to remind that these Mellin 
arguments are natural numbers. 

We start with introducing some notation. 

3.1.1. Functions L(z;aj Vk , . . . ,uJ Vl ijk, ■ ■ ■ > ji)- Choose k > 1; Vk, ■ ■ • ,vi G 
V, and nonnegative integers jk,---,ji] h is convenient to add jo = 0. In our 
applications, j a will satisfy the following restrictions: 

ja <m Va -1+ja-i- (3.2) 

Now put 

L(z;u> Vk , . . .,LO Vl ;j k , ■ ■ -,ji) := 
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r r p 1-Ki(n x ^ Vn k )z 



n™ vl+jo jl (m + n 2 ) m -2+n-j2 . . . ( ni + . . . + n ^m Vk +j k . 1 - Jk 



n 1 ,...,n k >l IL 1 

(3-3) 

Thanks to the presence of exponential terms in (3.2), this series absolutely converges 
for any z with Imz > and defines a holomorphic function in H. 

Notice that the enumeration of arguments of L is reversed in order to get a more 
natural enumeration of factors in the summands of (3.3). 

3.1.2. Numbers L(0; u> Vk , . . . ,co Vl ; j k , jk-i, ■ ■ • If we formally put z = 
in the expansion for 

(2niz) ~ 3k L(z;u Vk ,...,uJ Vl ;j k ,..., ji) 
we will get the formal series 

Cy 1 ,n 1 ■ ■ ■ c v k ,n k .\ 

m,_i>i n? vi+j0 ~ jl (ni + n 2 ) m ^+h-h . . . ( ni + . . . + nfe )m Vk+jk _ 1 - jk ' 
We have 

c Vl ,m ■ ■ ■ c Vk , nk = 0((mn 2 . . . n k ) c ) = 0((m + ■ ■ ■ + n k ) kC ) . (3.5) 
Assume that (3.2) holds. Then the general term of (3.4) is bounded by 

1 

mim + n 2 ) . . . (m + • • • + n fc _i)(ni + • • • + n k ) m ^ + ^-^- 1 - kC ' 

Hence (3.4) absolutely converges as long as 

m Vk +j k -i -jk>l + kC. 

Summarizing, we get three alternatives, describing the possible behavior of 

L(z; u> Vk , . . . ,u Vl ; j k , jk-i, ■ ■ ■ , ji) 

as z — > 0. We will later identify the respective limit as a component of the total 
Mellin transform of (ov): 

Case 1: j k = and m Vk + jk-i > 1 + kC. Then the limit exists, and equals to 
the "multiple Dirichlet series" 

L(0; 0J Vk ,..., w Vl ;0, jk-i, ■ ■ ■ , ji) = 
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V- C "l,"l • • • C Vk,n k /g g\ 

ni ,^ k >in? vl+jo ~ j \n 1 + n 2 ) m "> + 3*-»...(n 1 + ' 

(jo = jk = appear at the r.h.s. only for the uniformity of notation). 

Case U: jk > and m„ fe +jk-i —jk > 1 + kC. Then the limit exists, and vanishes 
thanks to the factor (2nizy k : 

L(0;u Vk , . . . ,u Vl ; j k , jk-i, ■ ■ ■ , ji) = (3.7) 

Case 3: jk > and m Vk + jk-i — jk < 1 + kC. In this case an additional study 
is needed. 

We can now formulate the first main result of this section. 

3.2. Theorem. For any k > 1, (v\, . . . , Vk) G V k , and Imz > we have 



ra n -lm„ 2 -l+]'i m Vk —l+jk-i 

( _i)Ef =1 (^-D y: e ••• e (- 1 ) 

.71=0 J2=0 jfe=0 



Jfc x 



m„ x -1 ! m„ 2 -1+ji !... m„ fe -l+j fc _i ! . 

x -r— : —. L(z; u v , . .., u Vl ;jk, . . . , ji) . 

Ji'jr- ■ ■ -Jk' 

(3.8) 

The proof requires an auxiliary construction. 

3.3. Auxiliary polynomials D^^i n ^ m (t). Choose now k > 1; Vk, ■ ■ ■ , v\ e 

V, and positive integers nfc, . . . , ni. It is convenient to agree that for k = the 
respective families are empty. 

Define inductively polynomials 
putting = 1, and 

D^ 1 '"'"** 1 (t) = (l + d t )- 1 (D%>-> nk m ( ni + --- + n fe t \ . t m, k+1 -^ (39) 
m vi ,...,m Vk+1 \ / v <-/ v m vl ,...,m Vk i ^ _|_ _ _ _ _|_ / / ■> \ j 
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where 

(l + a,)" 1 :=£(-l) fc # 

£;>0 

as a linear operator on polynomials. 
For example, 

m vi —1 

Dz 1 (t) = (-ir^- 1 (m vi -i)\ jr -^fr 1 - (s- 10 ) 

31=0 Jl ' 

In particular, D^fi (0) = (— l) m "i~ 1 (m„ 1 — 1)! Furthermore, 

m "l~ f_iyi „.?i«'i+to„ 2 -1 

^ 1 ^ 2 w = (-ir-- 1 (m„ 1 -i)! E L ~r( 1 +^)" 1 ( 1 ni+n2)J1 = 

m ' 1_1 (-l)*nf m \^ +Jl (-l)^ 



(-ir^^-^-i).^-!)! e e 



J1= o ii ! K+n 2 )^ ^ j 2 ! 

(3.11) 

In particular, 

Dx$*,<f» = (-ir".- 1+ "*--v«,-i)!K,-i)! E . i , l ( „/ +Ti ; )jl 

The general formula looks as follows: 
3.3.1. Proposition. We have for k > 1: 

m vi —l m V2 — l+ji m Vk —l+j k - 1 

^v:X fc (*) = (-i) Ef = i(m -«- 1) E E ••• E 

(m^ -1)!K -l+ji)!...(m^ fc -l+j fc _i)! 

• • • jfc! 

x 1 fjk 

n'^iri! + my*-** . . . (ni H + nfc-i)^- 2- ^- 1 (m H + re*.)**- 1 

(3.12) 
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Proof. We argue by induction on k. Assume that (3.12) holds for k and apply the 
operator at the right hand side of (3.9) to the right hand side of (3.12). Looking 
for brevity only at the last line of (3.11), we get: 

— (i+d t )-H mv ^- 1+jk 

n x n (m + n 2 ) J1 - j2 . . . (ni H + n k )^-i-3k ( ni H + n k+l )^ 

1 

= — — x 

n 1 31 (m + n 2 ) J1 - j2 . . . (ni H h n k )^-i-Jk ( ni _| 1_ n k +i) Jk 

x m ^ 1+Jk (-iy^(-ir^- 1+jk (m Vk+1 - 1 + h)\ ti ^ 

Combining this with (3.12) for k, we get (3.12) for k + 1. 

3.4. Proof of Theorem 3.2. By induction on k we will prove the following 
formula: 

(27Ti) m ^+- +m ^lU"v k ,...,0 Jvi ) = 
. . ,~ r „27ri(niH hn fc )« 

= E rn" Vl ' n ^-\: k ' nk ( - - ^ D^- ; ^ (2m(n 1+ . ■ .+n k ) z). 
r, „ >i n i ( n i + n 2) V1 . . . (ni + • • • + n fc ) m ^ 

ffcl , .... /it. ^ _L 



ni,...,n fe >] 

Combining it with (3.12) and (3.3) we will get (3.8). 
For k = 1 we check (3.13) directly: 



(3.13) 



^r-J^KJ = (2™)™- f Wl>1 (z0 = 

</ ioo 

= (27ri) m «i ^ c vuni / e 2 ™ 121 ^ 1 " 1 ^!- 
m=i ^ io ° 



Putting in the ni-th summand t = 2niniZi, we can rewrite this as 



n f' 
ni =1 i- 



eH m ^- l dt. (3.14) 



Since 



J e*P(t)cft = e*(l + a t ) _1 P(t) + const , (3.15) 
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this is equivalent to (3.13). 

The inductive step from k to k + 1 is similar: assuming (3.13)fc, we have 



(2ni) m ^ + --- +m ^lUu Vk+1 ,...,u Vl ) = 

/ j c «fe+i,n fe + i / 
n fc+ i=l At 



(2ttz)^ +1 V c Vk+unk+1 I e 2 — ^r fe+1_1 x 



27ri(niH hn fe ) 2fe 

r, «>i n i (m+n2)™" 1 ...(nH + ^fc) 1 fe 



= (27Tz) m ^ + l 



C Vl ,n 1 ■ ■ ■ Cv k ,n k Cv k + 1 ,n k + 1 



ni ,..^ +1 >i (ni + na)™-! . . . (m + • • • + n*)^ X 

"2 



x / ^ ni+ -"^^ 1 ^D^- i ":^ (2m( ni + ■ ■ ■ + n k ) z k ) z^ h+1 ' dz k . 

J ioo 

(3.16) 

Putting here t = luiiyix + ■ ■ ■ + n k+ \) z k , we can rewrite the last integral as 



1 1 



(2yrz) m ^+i (m + • • • + n fc+ i) m ^+i 

that is, 



^ e t D nu-,n k ( U 1 + --- + n k \ t m Vk+1 -l dt 



• ' 27ri(mH hn fc+1 )2: x 



(27ri) m "fe+i (mH hn fc ) m "fe+i 

x (l + ^)- 1 (D- 7 ;;X fc ( , Wl + '" + , Wfc * m " fc+1 ^ 



t=2iri(n 1 -\ \-n k + 1 )z 



ni H hn fc+ i 

Substituting this into (3.16), we finally obtain (3.13) and (3.12). 
3.5. The limit z — > 0. We can try to get an expression for 

A°oo (^fc > • • • > ) 

as a (linear combination of) multiple Dirichlet series, by formally putting z = in 
the r.h.s. of (3.8). However, we will find out that this cannot be done automatically 
for a certain range of values of (j k ,j k -i), namely, for j k > m Vk + jk-i — 1 — kC: 
cf. Case 3 at the end of subsection 3.1.2. 



24 



To solve this problem, we will have for the first time to assume that z 1 ~ mv uj v (z) 
are of cusp modular type, say, for the group To(N), or for any modular subgroup 
which is normalized by the involution z i— > gz, 

fO -1 

We can then apply Proposition 2.4 which we reproduce and slightly augment: 

3.5.1. Proposition. Assume that coy as above is a basis of a space of 1- forms 
invariant with respect to g^. Then 

4W = (gN^iojv)))- 1 ^^). (3.17) 

Replacing the coefficients of the formal series at the r.h.s of (3.17) by their (conver- 
gent) representations via multiple Dirichlet series (3.8), we get such representations 

3.5.2. Main application. If we fix a modular subgroup normalized by g^, this 
proposition becomes applicable to any family of cusp forms of the type uJi(z)z m ~ 1 , 
m > 1, where uJi{z) runs over a basis of the space of forms of a fixed weight 2r, and 
m runs over [1, 2r — 1] (cf. (2.4)). Moreover, we can mix different weights, that is, 
take a finite union of such families. 

Passing to a different basis of such a space, we may even assume that (oj v ) 
consists of eigenforms for g N : g* N (u = ±1, for all s E V. 

i 

Coefficients of are the series (3.3) at z = -j=. 

§4. Shuffle relations 
between multiple Dirichlet series 

4.1. Notation. In this section, we will consider (formal) multiple Dirichlet 
series of a special form generalizing expressions (3.4), and deduce bilinear relations 
between them generalizing the well known harmonic shuffle relations involving shuf- 
fles with repetitions. 

Each such series will depend on a set of coefficients data C and several complex or 
formal arguments Si. Here are precise definitions. Let k > 1 be a natural number. 

4.1.1. Definition. (%) Coefficients data C of depth k is a family of numbers 
Cn,m indexed by two pairs of integers satisfying j > i > 0, j < k, and n > m > 0. 
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(ii) The multiple Dirichlet series associated with C and arguments si,...,Sk is 



T I \ Uk>j>i>0 C Uj,Ui 

L c (s 1 ,...,s k ):= ^ ~ ~ (4.1) 



o=mo<""i <---<u k ez 



Ui U ■ ■ - U 



k 



4.1.2. Examples, (a) Assume that Cn,m = 1 if m > or i > and put 
^ =a { " ) . Then 



C n,0 



„(1)„(2) 

Sl '-' Sfc)= S < Iff... Iff ^ 

is an usual multiple Dirichlet series. 

(b) Define Cy 77, ctS 111 3.1, and choose vi,...,Vk G F as in 3.1.1. Construct the 
coefficients data C putting 

C n,m ^ := c Vj,n-mi (4-3) 

and Cnjm = 1 otherwise. Then Lc(m Vl + jo — ji, . . . , m Vk + jk—i — jk) becomes the 
formal series (3.4) if we redenote Uj = ri\ + • • • + rij. 

4.1.3. Shuffles and a composition of the coefficients data. Let C = (cnjm) 

and D = (dn,m) be two data of depths k and / respectively. A (k, I, p) -shuffle with 
repetitions is a pair of strictly increasing maps a = (a\, 02), 

a i: [0,fc]-[0,p], a 2 : [0,Z]-[0,p] 

satisfying the following conditions: 

Pi (0) = (7 2 (0) = 0, (7! ( [0, k] ) U (72 ( [0, I] ) = [0, p] . (4.4) 

It follows that max (k,l) < p < k + l. We will say that the a -multiplicity of j G [0,p] 
is one, if j ^ cr 1 ([0, fc]) fl CT2QO, /]). Otherwise the a-multiplicity of j G [0,p] is two. 
In particular, the cr-multiplicity of is two. 

Given such C, D and a = ((71,(72), we will define the third coefficients data 
E = (en,m) of depth p which we denote E = C* a D. Choose j, i with p > j > i > 0. 
We have the following set of mutually exclusive and exhausting alternatives (A) 1; 
(A) 2 , (B), (C). 
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(A) Assume that both j and i have multiplicities one. 

(A)i. Both j and % belong to the image of one and the same a a with a = 1 or 
a = 2. Then we put 

eg'S :=c& 1(jWl(l)) fora = l, 

and 

e^:=dS (]Kl(,)) forfl = 2. 

(A) 2. Assume that j, resp. i, belongs to the image of cr a , resp. a&, with a ^ b. 
Then we put 

p 0',i) ._ i 

(B) Assume that exactly one of j, i has multiplicity two. Then there exists only 
one value a = 1 or 2 such that j and i belong to the image of a a . We put then as 
in the case (A)i 

e { li :=&^ 1{l)) fora = l, 



and 



p (j,i) ._ /°2 1 U),°2 "(0) f _ 2 



(C) Assume that both % and j have multiplicities two. Then we put 

4.1.4. Shuffles and a composition of the arguments. Let s := (si, . . . , s^) 
and £ := (t\, . . . , t{) be arguments for the data C and D as above, and a a (/c, /, p)- 
shuffle as above. We define s + CT t := (ri, . . . , r p ) as follows. 

If z has multiplicity one and is covered by cxi, resp. 02, then 7^ := s^-i^, resp. 

If z has multiplicity two, then ri := s a -i^ +t a -~ 1 ( i y 
We can now state the main result of this section. 

4.2. Theorem. Let C , resp. D , be some coefficients data of depths k, resp. I, 
as above. Then we have 

L c (s)-L D (t) =J2 L c*„D(s+at) (4.5) 
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where the summation is taken over all (k, I, p) -shuffles with repetitions. 

Proof. Consider a term of the series (4.1) corresponding to (uo = 0, u\, u 2 , ■ ■ ■ , Uk) 
and a term of the series Lo(t) corresponding to, say, (wq = 0, w\, w 2 , ■ ■ ■ , wi). This 
pair of terms determines a unique (fc, Z,p)-shuffle (ui, a 2 ), where p is the cardinality 
of the union of sets 

{u 1 ,u 2 ,...,u k }U{w 1 ,w 2 ,...,w l } := {<?i,..., <? P }. 

Namely, we may and will assume that go = < qi < ■ ■ ■ < q p . Then = j if 

Ui = qj, and a 2 {i) = j if Wi = qj. 

Group together all pairwise products corresponding to one and the same shuffle, 
and denote the resulting sum by L a . 

The denominator of one such a product will obviously be q^ 1 . . . qp p where r = 
s + a t. Moreover, knowing such a denominator, we uniquely reconstruct the two 
terms from Lc(s) and Lo(t), from which it was produced, at least if s, t take generic 
values so that in the family {s a , tf,, s a + 1^} all terms are pairwise distinct. Finally, 
all possible sequences qo = < qi < ■ ■ ■ < q p will occur. 

To prove that L a = Lc , * ct d(s +a t), it remains to check that the numerator of 
such a product will be as predicted by (4.5), in other words, that 

p>j>i>0 k>j>i>0 l>j>i>0 

if en,m are defined as in 4.1.3. 

This is straightforward, although somewhat tedious. 

4.3. Concluding remarks. It would be interesting to describe some nontrivial 
spaces of Dirichlet series containing periods of cusp forms, closed with respect to 
the series shuffle relations, and consisting entirely of periods in the sense of [KonZa] . 

Regarding shuffle relations themselves, motivic philosophy predicts that they 
should be obtainable by standard manipulations with integrals. For the harmonic 
shuffle relations between multiple zeta values, A. Goncharov established this in the 
Ch. 9 of [G06] (for convergent integrals), and in 7.5 of [G06] elaborating the last 
page of [Go4] (for regularized integrals). Conversely, integral shuffle relations can 
be deduced from harmonic ones: see [Go5], Ch. 2. 



2<S 



§5. Iterated Eichler— Shimura and Hecke relations 

5.1. Eichler— Shimura relations for iterated integrals. In this subsection, 
we take for X the upper half-plane H and for (u>v) a family of 1-forms of cusp 
modular type (see 2.1.1) spanning a finite-dimensional linear space stable with 
respect to the modular transformations 

-1\ _fl -1 

1 J ' T \1 

5.1.1. Proposition. With these assumptions, we have 

<r.(J* (u v ))-J* (<j v ) = l. (5.1) 

r,VrK)) • t.(J£°(uv)) ■ JtM = 1. (5.2) 
Proof. From (1.15) we get 

which shows (5.1). 

Similarly, r transforms (0, zoo) into (zoo, 1), then to (1, 0), and according to (1.9) 

J Q 1 {uv)-JL{"v)-Jir{uv) = l (5.3) 

which is the same as (5.2). 

Notice that some care is needed in establishing (5.3): the geodesic triangle with 
vertices 0, zoo, 1 should be first replaced by a sequence of geodesic hexagons lying 
entirely in H and cutting the corners of the triangle, and then it must be checked 
that in the limit the hexagon relation replacing (5.3) tends to (5.3). This is routine 
for cusp modular 1-forms, cf. 2.1.3. 

We now pass to the relations involving Hecke operators. 

5.2. Hecke operators. In this subsection, V denotes a finite set, uj v = f v (z)dz 
a family of modular forms of weight two, and p v a family of primes, both indexed 
by v & V. Moreover, we assume that T Pv iv v = \ v co v , X v G C, where T Pv is the 
Hecke operator 

V 7 b=0 V 7 6=0 
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Put U := V \J V where V is another copy of the indexing set V, and for v' G V 
corresponding to v G V put uv := (p v )*(a;„). Let u>u be the family consisting of 
all ui v and uv- When we consider formal series of the type J^(uu) as in §1, we 
denote the variables corresponding to V, resp. V' , by A in resp. A v '. 

Denote by W the set of pairs w = (v, b) where v G V and b G [0, p v — 1] . Let uw 
be the family consisting of LO( v ,b) '■= h(p v , b)*(u> v ). When we consider formal series 
of the type J*(u>w), we denote the variables corresponding to w, by B w or B( v b y 

Define the following two continuous homomorphisms of rings of formal series: 

/ : C((Au)) - C((A V )) : l(A v ) := X V A V , 1{A V ,) := -A v , (5.5) 
r : C((B W )) ^ C((A V )) : r(B {Vjb) ) := A v . (5.6) 
5.3. Theorem. We have 

I(JLM) = r(JLM) • (5.7) 

Proof. We will check that 

I(JLM) = JL((^ V - u v ,)) (5.8) 

whereas 

r(J? QO (u w )) = J£o((E w («.»))) ^ 

6=0 

where at the right hand sides we consider both families as indexed by V. Since 
from (5.4) and the definitions above we obtain for each v G V 

XyUJy - U V > = ^(«,6)- (5.10) 

6 = 

this will prove the theorem. 
We have 

oo 

JLM = J2 Yl ^•••\4kr.,w t(1 ). (5.11) 

n=0 ( Ul ,...,u n )EU n 
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Consider one summand in (5.11). In the sequence (ui, ...,%) there are several, say 
< k < n, elements Vi G V and the remaining n — k elements Vj £ V . Application 
of / eliminates all primes in the subscripts of A Un . . . A Ul and produces a monomial 
in A v ; besides, it multiplies this monomial by (— l) n_fc Y[ Hence the coefficient 
at any monomial A Vn . . .A Vl in l^Jf^uu)) can be written as 

e (- i ) n_isi nMt^k,.,^)) (5.i2) 

sc[i,...n] \ies J 

where the operator its replaces Vj by v'j whenever j ^ S. 

On the other hand, the similar coefficient in J^^XyUy — uv)) is 

Iioo(^v„Vv n - U v > n , . . . , X Vl UJ Vl - Uy'J 

which obviously coincides with (5.12) because the iterated integrals are polylinear 
in ui. This proves (5.8). 

The check of (5.9) is similar. We have 

= E B (v n ,b n )--- B (v 1 ,b 1 ) I ioo( UJ (v n ,b n ),---^(u 1 ,b 1 ))- 

((v n ,b n ),...,(v 1 ,b 1 ))ew n 

Application of r produces a series in (A v ) whose coefficient at A Vn . . . A Vl equals 

(6i,...,b n ) 

This is the same as the respective coefficient at the r.h.s. of (5.9). 

§6. Differentials of the third kind 
and generalized associators 

6.1. Normalized horizontal sections. In this subsection, following [Dr2], 
we will define and study solutions of the differential equation (1.4), dJ z {Vt) = 
Q(z) J z {Vt) in the case when O = A v u v may have a logarithmic singularity at 
a point a so that it cannot be normalized by the condition J a (0) = 1 and in fact 
cannot be defined by the series (1.3). 
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We start with a local situation. Put 

r v>a := res a u v , R a := res a O = ^ r v , a A v . (6.1) 

V 

The normalized solution will depend on the choice of a local parameter t a at a, and 
a branch of logarithm log£ a . 

Let U be a disc around a uniformized by t a . Denote by logt a the branch of 
the logarithm in U which is real on lmt a = 0, Ket a > 0. Delete from U a cut 
from a to the boundary which does not intersect the latter interval and denote U' 
the remaining domain. Write t Ra for e Ralogta . It is a formal series in A v with 
coefficients which are holomorphic functions near a in U' . Assume that outside of 
a, all uj v are regular in U. 

6.1.1. Definition. A V a, -horizontal section J in U' is called normalized at a 
(with respect to a choice of t a ) if it is of the form J = K ■ t Ra , where K can be 
extended to a holomorphic section in some neighborhood of a in U which takes the 
value 1 at a. 

We will see that this definition produces a version of J%(Q). In fact, we get 
precisely if R a = so that t Ra = 1. 

6.1.2. Proposition. For any a andt a as above, there exists a unique normalized 
at a local section holomorphic in U' . 

Proof. In the course of proof, we will be considering only one point a, so we 
will omit it in the notation for brevity and write R, t, r v etc in place of former 

Ra: t ai r Vja . 

The equation Vn(if • t R ) = is equivalent to 

dK = tt'K + t~ 1 [R,K]dt (6.2) 

where 

Q' := O — R — = r v v v A v . 
t 

V 

We look for a solution to (6.2) of the form 

oo 

K = 1 + y fv 1 ,...,v n A Vl . . . A Vn 

n=l ( Vl ,...,v n )eV n 
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where f Vlj ... jVn must be holomorphic functions defined in some common neighbor- 
hood of a in U and vanishing at a. From (6.2) we find df v = v v so that the only 
choice is f v (z) := J z v v . Notice that all v v are regular in a common neighborhood 
of a uniformized by t a . We will see that in this neighborhood all other coefficients 
can be defined as well. 

If fv!,...,v„-i with required properties are defined for some n — 1 > 1 , we find 
from (6.2) ' 

dfv!,...,V n = h/ V 1 fv 2 ,---,V n +^ ( r v 1 fv 2 ,---,V n — fv2,...,V n -i r v n ) dt. 

By the inductive assumption, the r.h.s. is well defined and regular at a, so inte- 
grating it from a to z we get f Vl ,...,v n ( z )- 

6.2. Scattering operators. Let now a, b be two points where O may have 
logarithmic singularities. Choose t a , h as above, construct the neighborhoods U a , 
Ub and neighborhoods with deleted cuts U' a , U' h in which we have the holomorphic 
normalized horizontal sections J a , Jb- Now embed U' a , U' b into a connected simply 
connected domain W to which both J a and Jb can be analytically extended. Clearly, 
they are invertible elements of O x ((A v )} at almost all points x G W. Put 

J a h =J- x Jb. (6.3) 

As in the proof of Proposition 1.2, one sees that J\ G C((Ay)). Borrowing the 
physics terminology, we can call this transition element the scattering operator. 

I added twiddle in the notation in order to remind the reader that J h a , besides fi, 
depends on t a and tb as well, if at least one of the residues R a , Rb does not vanish. 
This dependence however is pretty mild. Let J' a = K'{t' a ) Ra be another horizontal 
section normalized with respect to some t' a . Denote by r a G C* the value of t' a /t a 
at a, and let t' a = T a ■ t a r a , T a (a) = 1. 

6.2.1. Proposition, (i) We have 

J' a = Ja-T*% K = K'-T R «. (6.4) 
(ii) Therefore, after replacing two uniformizing parameters t a ,tb by t' a ,t' b , we get 

j'a = T -R a ja T R K (6 _ 5) 
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Proof, (i) We have 

J' a = K'(t' a ) Ra = K' ■ (t a ) R "r R ° 

Since r Ra G C{{Ay)), and this section is invertible, K' -T Ra •{t a ) Ra is V^-horizontal 
as well, and since K' ■ T Ra at a equals 1, it is normalized, so K' ■ T Ra ■ (t a ) Ra = J a , 
which proves (6.4). 

(ii) Similarly, from J' b = J' a J' b a and (6.4) we get 

L — T T Ra J' a r~ Rb 

which together with (6.3) proves (6.5). 

6.3. Example: Drinfeld's associator. Let X = P^C), V = {0, 1}, 

1 dz 1 dz 

2ni z ' 1 2ni z — 1 

Then 

O = Aqujo + A\U)\ 

has poles at 0, 1, oo with residues Aq/2-ki, Ai/2ni, — (Ao+Ai)/2ni respectively. Put 
to = z, t\ = 1 — z. Then Jq in our notation is the Drinfeld associator <j>Kz(Ao, A\) 
from §2 of [Dr2]. 

6.4. Generalized associators. An essential feature of the the last example is 
that uj v are global logarithmic differentials on the compact Riemannian surface P 1 . 

Generally, let X be such a surface, and (a^) a finite set of iV points on it. The 
dimension of the space of global logarithmic differentials with poles in this set, 
Y^Cidlog fi, where q G C, /$ are meromorphic on X, is bounded by N — 1. It 
achieves the maximum value N — 1 iff the difference of any two points ai — cij is 
torsion in the divisor class group. I will call such a set (a^) logarithmic. The supply 
of logarithmic set depends on the genus of X. 

a) Genus zero. Any finite set of points is logarithmic. The respective iterated 
integrals in this case include multiple polylogarithms introduced in [Gol]. In fact, 
as Goncharov remarked, general iterated integrals can be reduced to multiple poly- 
logarithms. 

b) Genus one. In this case we can take any finite set of points of finite order, for 
example, the subgroup of all points of a given order M . 
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For a general subset of points, Goncharov found a Feynman integral presentation 
(in the sense of the last section of [Go2]) of the respective real periods. His formulas 
involve the generalized Kronecker-Eisenstein series and can be considered as an 
elliptic version of multiple zeta values. A similar problem is addressed in the recent 
work of A. Levin and G. Racinet. 

c) Higher genera. If genus of X is > 1, the order of such a logarithmic set is 
bounded. The most interesting explicitly known examples are modular curves and 
cusps on them, cf. [Drl] , [El] . 

Notice, that the initial Drinfeld's setting is modular as well: P 1 with three 
marked points "is" the modular curve r (4) \ H together with its cusps. 

According to Deligne and Elkik ( [El] ) , a set of points is logarithmic iff the mixed 
Hodge structure on Q) (where X° is the complement to ) is split that is, 

the direct sum of pure Hodge structures. 

6.4.1. Definition. Let X be a compact Riemannian surface, and (a^) a loga- 
rithmic set of points on it. Then any scattering operator of the form J a - is called 
a generalized associator. 

6.5. Relations between scattering operators. Three types of relations 
established for J*(u)y) in §1 are extended below to the case of the general scattering 
operators. 

6.5.1. Group-like property. We have 

A(J h a ) = J h a ®J h a (6.6) 

where A is defined in 1.4-1. 

To see this, notice that A(J ), by definition, is the series which is obtained from 
J a by replacing each A v with B v := A v ® 1 + 1<S)A V . Hence A(R a ) = R a (g> 1 + 1 ®R a 
is the residue of A(O) at a. Therefore A(J a ) is the normalized VA(n)~horizontal 

section in the ring of formal series in B Vl and A(J%) is the respective scattering 
operator. 

On the other hand, J a <S> J a satisfies the same equation d(J a ® J a) = (O <8> 1 + 1 <E> 
£i)(J a <S> J a ) and clearly is as well normalized at t a . Hence the passage from J a ® J a 
to Jb <8> Jb is governed by the same scattering operator, this time represented as the 
r.h.s. of (6.6). 

6.5.2. Cycle identities. Let 7 be a closed oriented contractible contour in 
U , inside which there are no singularities of O. Let a\, . . . , a n be points along this 
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contour (cyclically) ordered compatibly with orientation. Then 

J a a 1 J a al---J a a :- 1 J a a ; = ^- (6-7) 

Of course, in this statement we assume that at each point a\ one and the same 
local parameter t ai and one and the same branch of its logarithm is used for the 
definition of the two relevant normalized sections corresponding to the incoming 
and outcoming segments of the contour. Otherwise the relevant r R factors as in 
(6.5) must be inserted. 

6.5.3. Functoriality. Let g be an automorphism of X compatibly acting upon 
all the relevant objects, with the possible exception of parameters t a , and transform- 
ing the space spanned by u v into itself. Define g* as in 1.4- Then 

J 9 g b a = T^g*(J b a )r b R ». (6.8) 

where the t r factors account for the passage from t a , resp. t^, to g*(t a ), resp. 

9*(ta). 

The proof is essentially the same as in 1.4 and 1.5. 

6.6. Example: Drinfeld associator revisited. If we treat Drinfeld's setup 
as the ro(4)-modular curve, lift it to H and apply to the respective family of 
scattering operators the Eichler-Shimura relations (5.1), (5.2) (following from the 
cycle identities of lengths 2 and 3, and functoriality), we will get the duality and 
the hexagonal relations. Of course, this is how they were deduced in the first place, 
with cr, r pushed down to P 1 rather than everything else lifted to H. 

It seems very probable that the somewhat mysterious relationship between the 
double logarithms at roots of unity and the modular complex discovered in [Go3] 
can be explained in the same way. 
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